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Abstract 

Periodic orbit theory is an effective tool for the analysis of classical and 
quantum chaotic systems. In this paper we extend this approach to stochas- 
tic systems, in particular to mappings with additive noise. The theory is 
cast in the standard field theoretic formalism, and weak noise perturbation 
theory written in terms of Feynman diagrams. The result is a stochastic 
analog of the next-to-leading Ti corrections to the Gutzwiller trace formula, 
with long time averages calculated from periodic orbits of the deterministic 
system. The perturbative corrections are computed analytically and tested 
numerically on a simple 1-dimensional system. 



1 Introduction 



Noise plays important role in a variety of physical contexts. Robustness to noise is 
of interest for any system since there is always some small length scale at which the 
dynamics is affected by thermal or quantum fluctuations or unobserved degrees 
of freedom. For example, the interplay of deterministic dynamics and magnetic 
diffusivity is subject of great interest in the dynamo problem, where the effect of 
magnetic field diffusion on the steady fast kinematic dynamo rates is discussed 
in |I| within the periodic orbit theory formulation of refs. |2|, ||, f|. 
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The noise tends to regularize the theory, replacing the deterministic delta 
function evolution operators by smooth distributions. While in this paper we are 
interested in effects of weak but finite noise, the a — » limit is also important as 
a tool for identifying the natural measure |5|, |(| [7j for deterministic flows. The 
noise regularization might in addition cure some of the ills of intermittent systems 
which are plagued by power-law convergences arising from terms like |A — 
in the limit A — > 1. 



We have cast the theory in the standard field theoretic language [||, in the 
spirit of approaches such as the Martin-Siggia-Rose |J formalism, the Parisi- 
Wu ju|] stochastic quantization, and the Feigenbaum and Hasslacher Q] study 
of noise renormalization in period doubling. This perturbation theory has the 
same structure as the % corrections to the semiclassical Gutzwiller trace formu- 



las JT^] computed by Gaspard and Alonso [13, 14, |15|| , and the trace formulas for 



continuous stochastic flows and for the Ti corrections formulated by Vattay [1 

Though it is clear from the literature on stochastic path integrals that some 
kind of Feynman diagrams apply, the present work seems to be one of the few that 
actually compute the weak noise corrections for a concrete dynamical system, 
although in some cases the leading correction may be obtained directly from 
the perturbed eigenfunction [[17], 18 1. The form of the perturbative expansions 



of sect. |4| is reminiscent of perturbative calculations of field thery, but in some 
aspects the calculations undertaken here are relatively more difficult. The main 
difference is that there is no translational invariance along the chain, so unlike 
the case of usual field theory, the propagator is not diagonalized by a Fourier 
transform. We do our computations in configuration coordinates. Unlike the most 
field-theoretic literature, we are neither "quantizing" around a trivial vacuum, nor 
a countable infinity of stable soliton saddles, but around an infinity of nontrivial 
unstable hyperbolic saddles. 

Two aspects of our results are a priori far from obvious: (a) that the structure 
of the periodic orbit theory should survive introduction of noise, and (b) a more 
subtle and surprising result, repeats of prime cycles can be resummed and theory 
reduced to the dynamical zeta functions and spectral determinants of the same 
form as the for the deterministic systems. 

Having constructed the perturbation expansion in sect. |], in sect. |5| we con- 
front the theory with a numerical determination of eigenfunctions and eigenvalues, 
and verify the correctness of our perturbation expansion to the same numerical 
accuracy. A variety of flow models with noise are simpler to study in nonper- 
turbative (large a) limits; numerical eigenfunctions do not depend on the weak 
noise assumptions, and in fact require the noise to be larger than the effective 
discretization length of the basis. 
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2 Stochastic evolution operator 



The periodic orbit theory allows us to calculate long time averages in a chaotic 
system as expansions in terms of the periodic orbits (cycles) of the system. The 
simplest example is provided by the Perron-Frobenius operator 



£>p{y)= / dx5{y - f{x))p{x) 



for a deterministic map f(x) which maps a density distribution p(x) forward in 
time. The periodic orbit theory relates the spectrum of this operator and its 
weighted evolution operator generalizations to the periodic orbits via trace for- 



mulas, dynamical zeta functions and spectral determinants [19, |20|. Our purpose 
here is to develop the parallel theory for stochastic dynamics, given by the discrete 
Langevin equation [^l], ^] 

Xn+l = f(x n ) + er£ n , (1) 

where the £ n are independent normalized Gaussian random variables. 

We shall treat a chaotic system with such Gaussian weak external noise by 
replacing the the deterministic evolution 5-function kernel by C a , the Fokker- 
Planck kernel corresponding to (||), a sharply peaked noise distribution function 

C„ = 5 a (y-f(x)), (2) 

where 5 a is the Gaussian kernel 



The method can be applied to smooth distributions other than the Gaussian one 
in the same manner. 

We shall evaluate the trace formulas by steepest descent methods, and obtain 
the noisy traces (traces of C a ) and determinants in terms of the cycles of the 
deterministic system. The theory is then tested numerically on one-dimensional 
maps, but we expect the generalization to higher dimensions to be of the same 
structure as the formulas derived here. 
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3 Stochastic trace formula, steepest descent approxi- 
mation 



We start by calculating the trace of the nth iterate of the stochastic evolution 
operator L a for a one-dimensional analytic map fix) with additive Gaussian noise 
a. This trace is an n-dimensional integral on n points along a discrete periodic 
chain, so x becomes an n-vector x a with indices a,b, . . . ranging from to n— 1 
in a cyclic fashion 

n—l t_ 

x n = x , [dx] = Y[ a . (4) 

a=0 V2TT(T 2 



As we are dealing with a path integral on a finite discrete chain, we find it 
convenient to rewrite the exponent in matrix notation 

tr£- = J [dx] e-lA-WM 2 /^ ; Kb = 6ab+1 7 (5) 

where x and f(x) are column vectors with components x a and f(x a ) respectively, 
and h is the left cyclic shift or hopping matrix satisfying h n = 1, h~ l = h T . Un- 
less stated otherwise, we shall assume the repeated index summation convention 
throughout, and that the Kronecker 5 function is the periodic one, defined by 



I it — _i_ 

For sufficiently short chains, (Q) is an integral that conceivably lends itself 
to numerical evaluation [23], although clearly not in the long time n — > oo limit. 
However, if the noise is weak, the path integral (||) is dominated by periodic 
deterministic trajectories. Assuming that the periodic points of given finite period 
n are isolated and the trajectory broadening a sufficiently small so that they 
remain clearly separated, the dominant contributions come from neighborhoods 
of periodic points; in the saddlepoint approximation the trace (||) is given by 

trC^ £ e w °, (7) 
x c eFix/" 
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where the sum goes over all periodic points x c = x c+n of period n, f n (x c ) = x c . 
The contribution of the x c neighborhood is obtained by shifting the origin of 
integration to 

where from now on x a refers to the position of the a-th periodic point, and 
expanding / in Taylor series around each of the periodic points in the orbit of x c . 

The contribution of the neighborhood of the periodic point x c is given by 

e W c = J [d(j)]e -(A-^V> W j 2 /2^ 

= |detA| J '[dip] e£i tr (AV"(0))V* 2 /2* 2 ( 8 ) 
where the propagator and interaction terms are collected in 

oo f h m + 1 



a m=2 



We find it convenient to also introduce a bidirectional propagator C = AA T 
for reasons that will become apparent below. In the second line of ^ we have 
changed coordinates, 

<p = A" V - V'W , (10) 
and used the matrix identity In detM = tr In M on the Jacobian 

I - - detA e" tr "n(i-AV") (in 

det (A-i - V") ~ det (1 - AV") ~ ' 1 j 

The functional dependence of = (f>(ip) is recovered by iterating fllOD 

4> a = A ab <p b + A ab V b '{4>) . (12) 

The above manipulations are standard |0] and often used in the "stochastic 
quantization" literature [10, where they are artfully employed to promote 



identities such as detM j 'detM = 1 to supersymmetric field theories. Such sym- 
metries do not seem to simplify the calculation at hand. 
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The saddlepoint expansion is most conveniently evaluated in terms of Feyn- 
mann diagrams, which we now introduce. The interaction terms in V and its 
derivatives can be represented in terms of the vertices 



and the propagators as directed lines 
A a6 = — «— 

C a b = * • * 

The first two derivatives of V may be written 




with the cross-hatched circle as V and the diagonally filled circle as <f>. The 
relation between the fields (|T^ ) becomes 




with the small open circle as (p. This recursively generates all tree diagrams 
ending in ip, which the stochastic averaging of sect. ^| will tie into loop corrections. 

As the sum is cyclic, e Wc is the same for all periodic points in a given cycle, 
independent of the choice of the starting point x c . 

In the saddlepoint approximation we assume that the map is analytic and 
the extrema f n are isolated. For the leading a 2 correction that we shall evaluate 
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here we need derivatives of / up to the third. A map with non-analytic points or 
marginal stability would lead to additional diffraction corrections that we shall 
not consider here. 



From the second path integral representation in @ it follows that A can 
be interpreted as the "free" propagator. As A will play a central role in what 
follows, we write its inverse in its full [nxn] matrix form: 



A" 1 = h- 1 - f 



-fx 1 



-/a 1 



V i 



\ 



-f'nJ 



(16) 



where f is a diagonal matrix with elements f' a = f'(x a ) a shorthand notation for 
stability of the map at the periodic point x a . The determinant of A is 



n-l 



det A 



A c -1 



Ac=]Jf\xa), 



(17) 



a=0 



with A c the stability of the n cycle going through the periodic point x c . We shall 
assume that we are dealing with a chaotic dynamical system, and that all cycles 
are unstable, |A C | > 1. 

The formula for propagator itself is obtained by inverting (|l^) and using 
relation (hi') n = A c , (due to the periodicity of the chain): 



A 



1 - f'- 1 /^ 1 



^(f'^vyf- 1 

fc=0 



1 



n-l 



k=0 



(18) 



In the full matrix form, the propagator is given by 



/ fl—fn-1 


fl---fn-\ fz—fn—\ 


1 \ 


1 


/2---/o /3/4—/0 


... / 0# 

/0/1 


A 


1 f-3---fofl 


/1/2 




/0/1/2 


/1/2/3 


f2f3 /3 




V /l---/n-2 


/2"-/n-2 


1 /o-../n-2/ 
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or, more compactly, 



A 



-1 ^ j, . . -I 



ah 



- J] / (xd) , A aia „ 1 = r — (20) 



where cf increases cyclically through the range 6 + 1 to a — 1; for example, if 
a = 0, a — 1 = n — 1. We note that A is invertible only for cycles which 
are not marginal, |A C | ^ 1. The |A C | = 1 case we would require going beyond 
the Gaussian saddlepoints studied here, and typically to the Airy-function type 
stationary points p5|. 



4 Weak noise perturbation expansion 

The saddlepoint approximation (||) is a discrete path integral on periodic chain of 
n points which we shall evaluate by standard field-theoretic methods. Separating 
the quadratic terms we obtain 

e w c = _i_ J [M e -s ^-s lM (21) 

where 

S (<p) = <p 2 /2a 2 , S^) = -E fc tr " (22) 

k=l 

The terms collected in Si(ip), linear or higher in (p, are the interaction vertices. 
Next introduce a source term J a and define a partition function 

e W c (J) = 1 f [ d¥ ,] e -So((p)-S/((p)+JalP« 

|A C 1| J 



A, -I 



1 



-Sdh) e ^ J \ (23) 



|A C - 1| 



Here we have used standard formulas for Gaussian integrals together with the 
normalization (||). In our diagrammatic notation this is 



e W c [J] = _^ e {( ^ + 2 T^f—lJ >• " (21) 



8 



Expanding 



CT 2 72 a 2 a A 

e- J =l + —J a Ja + ^JaJaJbJb + ..., (25) 
z o 



operating on this series with exp{— Sj(jj)}, 



d d 2 j t 

g — JdJd 



dJ a dJ b 



= cr 2 5 a b ■ 

J=0 



collecting terms of the same order in a 2 , and setting J a to zero yields the pertur- 
bation expansion 

oo 

W c = -ln\A c -l\ + Y,W c ,2k(T 2k ■ (26) 
fc=i 

In field-theoretic calculations the W c q term is usually an overall volume term 
that drops out in the expectation value computations. In contrast, here the 
W C) q = — In | A c — 1| term is the classical weight of the cycle which plays the key 
role both in the classical and stochastic trace formulas. 

In diagrammatic language, we join all possible pairs of ip vertices, each one 
giving a a 2 C propagator. Thus the first diagram in (H) is expanded (||jl|) to 



and then the <p vertices joined to form two diagrams contributing at order a 2 . The 
full noise corrections of order a 2 are given by all connected two-loop diagrams: 





W, 



c,2 



lOOiOOKX)^ (27> 



Each diagram has a two-fold symmetry, hence all combinatorial weights equal 
1/2. Before writing down the final expression, we note that several sub-diagrams 
may be simplified using (p0[). These are (no sum on a, b) 



a a _ A c 1 + (A c - l)5 ab 
ab ba ~ (Ac - I) 2 fj b 
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Adding the terms we obtain the a 1 contribution to the trace: 



2 A c -1 



1 



A, 



c 




aa 



+ 



Ag + A c 
(Ac - I) 2 



E Ja H fi 
pi j,i ^ab 

ab Ja Jb 



A c -1 



A c 



J2^fb^C bb . (29) 



ab Ja 



For an alternative approach to evaluating multiple derivatives, see appendix |A[ 



4.1 Repeats of prime cycles 

In the deterministic case repeats of periodic orbits can be summed up, and spec- 
tral determinants and dynamical zeta functions written in terms of prime cycles 
rather than periodic points. In order to accomplish this for the stochastic case, 
we need to compute the trace for repeats of periodic orbits. 

For r repeats of a prime cycle p we have n = n p r, A c = A p , where A p is 
the stability of the prime cycle p. Each index a = 0, . . . , n—1 is decomposed as 
a = a + an p + a with d = 0, . . . ,n p — 1 and a = 0, . . . , r— 1. d is an arbitrary 
starting point on the orbit which may be chosen independently for each index. / 
and its derivatives depend only on d. 



where a — bis the number of full repeats of p contained in nS=6+i' ^ ms ^ s achieved 
by setting d = b + 1. The sums over a and b are performed, leading to 



which is just a factor depending on r and A p multiplied by the sum for the single 
repeat of the prime cycle. 

When the calculations are carried out for both of the other sums, some rather 
unenlightening algebra leads to exactly the same prefactor; we discuss this rather 
remarkable point and its generalization to higher orders in detail in the sequel 



The first sum in (]29| ) is 
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paper |^] . Combining all three terms of (|29| ) leads to an expression for the trace 
in terms of cycles: 



l-zC„ ^ ^ 



~ z n P r J a 2 at (A r + j) | 
n P2^—r - eX Pi — W P2— ~ r ) ( 30 ) 

P 



n=l x c eFixf 

= E^EjArrii^iT^ (A--1) 2 

up to order cr 2 , where 

f" f" a 1 /" 

+ E -jr-jr C ab ~ A P 1 1 E ~~p~fb A a bC bb (31) 

ab Jo. h A P + 1 ab Ja 

contains all the dependence on the higher derivatives along the prime cycle p, 
with no dependence on the repetition number r. To put it another way, if p is a 
cycle, not necessarily prime, then 

Wpr 2 = rw p ,2 . (32) 

Next, using the identity 

00 



(1-x) 3 



k=0 



we rewrite the trace formula in a form in which repeats are resummed over by 
expanding the exponential in (|30| ) to order cr 2 , forming the sum over k, and 
putting the result back in an exponential: 



y n„r 2 



1 ZLj o V fc=0r=l l^pl^P 



00 



E^Et^T' (33) 



where t p j~ is the k-th. local eigenvalue 
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This is the stochastic equivalent of the Gutzwiller trace formula for the semiclas- 
sical case IT! 



We sum over r as usual [20] to obtain from (33) a Selberg type product for 
the noisy spectral determinant 



oo 

det(l - zC a ) = J] J] (1 - tp, k ) , (34) 
v k=0 

valid to order a 2 . 

We observe a crossover effect, since for higher order eigenvalues (large k), 
eventually the argument of the exponential becomes of order one, and further 
noise corrections are required. This is as it should be: the higher order eigen- 
functions have more detailed structure, are more quickly smeared by the noise, 
and should decay faster. 



4.2 Fixed point 

The contribution from a fixed point (cycle of length one) is particularly simple, 
as all the sums and products collapse to a single term, and / = A. We obtain 



£_)__£ 

A(A + 1) ] " I A - 1 J A - 1 



^ = T7rr7T 3 r^r " t^-t > (35) 



If the map contains only a single isolated unstable fixed point, we thus have 
an expression for the eigenvalues, 



1 1 



2 



z k {a) |A|A fc 

valid to order a 2 . Note that depending on the sign of w Pt 2 , small amounts of noise 
can either enhance or inhibit escape from the fixed point. Higher order terms for 
a fixed point are given in ref. [26|. 



5 Numerical tests 



To test the above expressions for the trace, we have computed the required deriva- 
tives for the 23 prime cycles up to length n = 6 for the quartic map 

/( r) 20 j _- Q-^Vj . (36) 
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The choice of the map is motivated by requiring that the system be simple (one- 
dimensional in this case), with non-trivial /", /"' (hence quartic), with complete 
binary dynamics (hence a nice repeller), and no diffraction and nonhyperbolic 
regions in the immediate vicinity of the repeller (where the Gaussian saddle points 
would be insufficient). 



5.1 Evaluation of the determinant 



In this we follow the approach to computing escape rates originally introduced by 
Kadanoff and Tang p7|]. The topological length truncated cycle expansions p0| 
of spectral determinant (^J) are obtained by writing the trace and determinant 
as power series expansions in z and a, 



zL„ 



tr— -z- = $> n (a^ + <r 2 c n , 2 ) (37) 

a n=l 

N 

det(l-z£ a ) = l-^z«( Cni0 + (J 2 c n , 2 ) (38) 



n=l 



Here, the C coefficients come from (j30|), and the c coefficients are obtained by 
equating coefficients in 

zC d 

det(l - zC a ) tr a — = -z— det(l - zC a ) (39) 



following from the identity In detM = tr In M. The solution is found recursively 
as 

Cn,0 = — [Cn,0 — Cn-l,0Cl,0 — • • • — C'l^Cn-l.o] (40) 

n 

C n ,2 = — [C n ,2 — (Cn-l,0Cl,2 + Cn-i^Cifl) — 

... - (Ci j0 C n _i i2 + C lt 2Cn-l,o)} (41) 



From the c n fi coefficients we construct the deterministic Fredholm determinant, 
from which the deterministic eigenvalue v§ is found using Newton's method on 
the characteristic equation for C a at a = 0: 



N 

E 

n=l 



^0 l °nfl 



(42) 
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n 


^0 


^0,2 


1 


0.308 


0.42 


2 


0.37140 


1.422 


3 


0.3711096 


1.43555 


4 


0.371110995255 


1.435811262 


5 


0.371110995234863 


1.43581124819737 


6 


0.371110995234863 


1.43581124819749 



Table 1: Significant digits of the leading deterministic eigenvalue and its a 2 co- 
efficient, calculated from the spectral determinant as function of the cycle trunca- 
tion length n. Note the super exponential convergence of both vq and 1^0,2 (n = 6 
result is limited by the machine precision). 



The <7 2 correction to the eigenvalue is found from the a 2 terms in the characteristic 
equation, and comes to 



EN ,,—n„ 
W=1 l/ Q C ra , 2 

Y^TV -n-l 



"2 - _ w _i — (43) 



The leading eigenvalue uq for the deterministic (noiseless) map and the coeffi- 
cient of the a 2 correction 1/0,2 > shown in table |], demonstrate the superexponen- 
tial convergence with n of spectral determinant, as expected for nice hyperbolic 
dynamical systems. 

The escape rate 7 of the repeller is calculated directly from the eigenvalue, 
7(0") = — In v(a) 
and hence 

70 = - In z/ 

72 = 

^0 

We have also directly tested the repeat formula (32) for our cycle set. 



5.2 Discretized eigenfunction 

In order to check the perturbative calculation we evaluate the eigenvalue numer- 
ically as a function of a by treating the evolution operator as a matrix acting on 
a discretized eigenfunction. That is, we approximate C by a matrix C ycXc where 
(x c ,y c ) is the center of a square in the x-y plane of small but finite size e (with 
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upper limits (xt, yt) and lower limits (cc&, Vb))- The matrix element is obtained by 
assuming the distribution p is constant across this small square: 

rx t ryt 

£y c x c = e / / 8 a (y - f(x))dxdy . 
Jxb Jyb 

The integral may be approximated using just a few evaluations of the kernel with 
errors of a similar order to those due to the variation in p. For example 

£-Vcx c = 2g(-^** + ^^tc + L t b + 2L ct + 4L CC + 2L c b + + 2L;, C + L^) 

with L tc = S a (yt — f(x c )), etc. requires only four evaluations per square since 
the boundary points belong to more than one square. For the map at hand, the 
discretized evolution operator leads to six digit accuracy in the escape rate for 
values of a as low as 3e. 

When a is small the matrix is very sparse, a fact which can be used to speed 
up the calculation. The leading eigenvalue is obtained by repeatedly evolving 
and rescaling an arbitrary smooth initial distribution, which then approaches the 
leading eigenfunction. In our case the map is expanding on the neighborhood of 
the asymptotic repeller, so the eigenfunctions are smooth, and the discretization 
procedure is stable. The numerical eigenfunctions for two values of a is shown 
in fig. |2[ The eigenfunction peaks at the critical point of the map, but this 
has no detectable effect on the eigenvalue, as subsequent iterations send the 
points around the peak towards — oo, away from the repeller x £ [0, 1]. As the 
flow is conserved only on the infinite interval x £ (— 00,00), no normalizable 
eigenfunction exists. On any finite interval, however, the escape rate is non- 
zero. If no point outside the interval can return (ignoring the exponentially 
small tunneling probabilities), that is, if the interval encloses the repeller and 
the neighboring region around it determined by the magnitude of the noise, the 
escape rate and the leading eigenfunction (up to normalization) are independent 
of the interval chosen. 

Subtracting the perturbative analytic terms from the numerically computed 
^o(cr), 

vq(<t) - 0.371110995234863 - 1.43581124819749a 2 ~ 0(a 4 ) (44) 

we compare with the numerically computed vq((t) in fig. [I], and estimate the next 
term to be approximately 38<r 4 . 
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difference 

10- 2 
10 " 3 
10 " 4 
10 " 5 

10 " 6 

0.01 0.02 0.03 0.05 0.07 0.1 

noise strength o 

Figure 1: The deviation of the analytic estimate ( |44|) from the numerically 
computed vq(g) for a range of values of the noise strength a (points), vs. the 
conjectured remainder, 38<r 4 (solid line). For small a the errors are due to the 
finiteness of the grid, and for larger a the deviation is due to the neglected higher 
order contributions. 



amplitude of 
eigenfunction 

10 
8 
6 
4 
2 


-0.5 0.5 1 1.5 2 

space coordinate 

Figure 2: The eigenfunction for a = 0.01 (sharp peak) and o = 0.1 (smoother). 
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6 Summary and outlook 



We have formulated weak noise perturbation theory for noisy maps in terms 
of periodic orbits of the deterministic system, expanding to order a 2 explicitly, 
resummed repeats of prime cycles, and tested the results numerically. From here, 
there are many possible generalizations and future directions. 



In the sequel paper [26] we shall recast the remarkable resummation of repeats 
(sect. 4J) in a more general framework applicable to all orders of the expansion. 
It seems from the numerical section that the coefficients of powers of a are growing 
very rapidly. The expansion is expected to be asymptotic. 

Evaluation of expectation values [20] on a stochastic flow requires replacing 
the Perron-Frobenius operator @) by the generalized evolution operator 

C <T = S a {y-f[x))eP< x K (45) 

The same general perturbation theory applies, but now an observable a{x) con- 
tributes an extra set of interaction vertices to Si(<ft) in J2l|). Similarly, the ad- 
dition of more dimensions and/or non-Gaussian weak noise can be treated by 
modifying the propagators and adding new vertices. 

While for deterministic flows it is appropriate to replace a flow by a return 
map on a Poincare section of the flow, it is not clear that this is appropriate 
for stochastic flows; a noise that is "white" on the Langevin equation level is 
"colored" when integrated to a Poincare section return, and it might have memory 
of the trajectory that a noisy iterated mapping cannot mimic. 

The noise in general has a different structure than the deterministic equations 
of motion; it typically breaks whatever symmetries the classical flow might have, 
unless clever precautions are taken to ensure that the noise respects the symme- 



try [28]. This situation is familiar from Quantum Mechanics, where quantization 



and canonical transformations do not commute. 

Our saddlepoint approximation to the spectrum of the exact evolution oper- 
ator receives perturbative contributions from all cycles, no matter how long. 
However, the noise causes the physical system to effectively lose memory at a 
rate depending on the region of phase space, so it might be possible to obtain 
accurate averages by replacing the evolution operator by effective finite memory, 
finite Markov partition transfer matrices. 

Such studies might enable us to understand the range of applicability of the 
"semi-classical" theory in greater detail than for the single cut-off time proposed 



in the case of semiclassical quantization by Berry and Keating [29]. Different 
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regions of phase space are dominated by different time scales, and the program 
of periodic orbit theory allows us to use the dynamics itself, encoded in the 
properties of cycles, to determine at what point classical behavior is modified by 
semiclassical or noise corrections. 



As in the semiclassical case, the saddlepoint approximation causes the multi- 
plicative structure of the evolution operators to be lost, and one might consider 
extended formulation of ref. [3C] to improve the analyticity of the spectral det- 
erminants. Finally, non-analytic points in the dynamics will lead to diffraction 
effects which are of different orders in <r, for example, the escape rate of the map 
4x(l — x) which has a quadratic maximum at the boundary of the deterministic 
repeller is of order 
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A Appendix: Recursive evaluation of derivatives 



The derivatives of x n = f n (x) 

i dx n n d x n in d x n (ac\ 

Xn ~lLx~ 1 Xn ~ lh? ' Xn ~ Ihf ' ' ' ' ' 1 ' 

with initial values 

dxo cPxq d 3 XQ 

xo = x, — = 1 , —y = , — r = , . . . , (47) 
ax ax z ax-* 

can be computed recursively by 

n 

x 'n+l = f'( x n) x n = Y[ f{Xk) 
k=0 

x n+l = / ( x n) ( x n) + / ( x n) x n = X n+ i f'(x k ) 
i3 



x n+l — / (^n) ( x n) '^f i x n) x n x n f i x n) x n 

, ^ f'"(x k ) , „ 2 , „ , X - f"( X j)f"( X k) i i 

x n+l = / ( x n) { x n) + 6/ (x n )x n (x n ) 

+4f"(x n Kx' n + 3/"(x n )«) 2 + f{x n )x'Z , (48) 



x' n+ i in the above has form of a propagator, f"(xk)/f'(, x k) 3- vertex, x'" +1 gets 
contribution from a 4-vertex diagram plus three 4-leg diagrams with two 3- 
vertices, etc.. In another words, this iteration of the "transport equations" gen- 
erates the Feynman diagram expansion. 
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